Finite-size effects on the dynamics of the zero-range process 
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We study finite-size effects on the dynamics of a one-dimensional zero-range process which shows 
a phase transition from a low-density disordered phase to a high-density condensed phase. The 
current fluctuations in the steady state show striking differences in the two phases. In the disordered 
phase, the variance of the integrated current shows damped oscillations in time due to the motion 
of fluctuations around the ring as a dissipating kinematic wave. In the condensed phase, this 
wave cannot propagate through the condensate, and the dynamics is dominated by the long-time 
relocation of the condensate from site to site. 
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The dynamics of fluctuations in simple noncquilibrium 
steady states of interacting particle systems has been 
studied extensively in recent years, and a fairly good un- 
derstanding of the physical processes involved has been 
achieved in infinite systems [1[ . Studies of these fluctu- 
ations in a finite system show a strong imprint of the 
nonequilibrium character, which combines with size ef- 
fects to bring in interesting dynamical phenomena 
Many of these systems often exhibit a nonequilibrium 
phase transition in the steady state, from a disordered 
to an ordered phase, as one tunes an external parame- 
ter such as the density. This raises the natural question: 
How is the dynamics of fluctuations in the steady state 
affected as the system passes through a nonequilibrium 
phase transition? In this paper, we address this issue 
within the ambit of a paradigmatic model, the zero-range 
process (ZRP) [1,0]) and show that there are strong dif- 
ferences in the dynamical properties arising from very 
different physical processes in the two phases. 

The ZRP involves biased hopping of particles on a pe- 
riodic lattice with a rate that depends on the occupancy 
at the departure site. At long times, the system reaches 
a nonequilibrium steady state. For certain classes of the 
hop rates, as the particle density crosses a critical thresh- 
old, there is a continuous phase transition from a disor- 
dered phase with uniform average density to a condensed 
phase where a finite fraction of particles (the "conden- 
sate" ) accumulates on a single site [4^] . The ZRP applies 
to a wide variety of physical systems ranging from traf- 
fic flow [H'] to shaken granular gases f6|; in addition, it 
was invoked to provide a criterion for phase separation 
in one-dimensional driven systems 0. 

Recent work on the ZRP has dealt with the relaxation 
of an initial homogeneous density distribution toward the 
condensed phase [8;]. By contrast, here we are interested 
in the dynamics of density fluctuations in the steady state 
in both the disordered and the condensed phases. We 
explore the dynamics by monitoring the variance of the 
integrated particle current in Monte Carlo simulations 
and supplement our flndings by analyzing the survival 
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FIG. 1: (Color online) (a) Schematic plot of the integrated 
current fluctuations C{t) across a bond in the steady state 
of the ZRP as a function of time t in the disordered phase 
(p < pc), at the critical point {p = pc), and in the condensed 
phase (p > Pc). (b) Schematic plot of C{t)/t as a function of 
t in all the phases. 



probability distribution of the largest mass in the system. 
The relevant time scales in the behavior of fluctuations in 
the two phases and the physical effects underlying them 
are summarized below. 

In Fig. [2 we show schematically the behavior of the 
variance of the integrated current as a function of time in 
different phases. In all the phases, at short times, the in- 
tegrated current is Poisson distributed, implying that the 
variance grows linearly in time, a behavior which holds 
for all times in an infinite system. In a finite system, in 
the disordered phase and at criticality, the variance shows 
oscillations at times proportional to the system size L. 
This results from kinematic waves transporting density 
fluctuations around the system with a well-defined speed. 
At longer times ('^ L^^"^), the wave decays, and then the 
variance increases linearly with time with a small slope 
that decreases with increasing system size. In the con- 
densed phase, however, the kinematic wave cannot pass 
through the condensate; thus, fluctuations do not circu- 
late around. The initial linear behavior continues until, 
after a characteristic time which grows as a power of the 
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system size, the condensate relocates itself. This results 
in the variance showing a linear rise in time with a much 
larger slope than at early times. Subsequently, after the 
condensate has relocated to another site, the slope of the 
linear rise slowly approaches a size-dependent constant. 

We now turn to a derivation of these properties. The 
ZRP involves N particles of unit mass on a ring of size L 
with arbitrary occupancy allowed at any site. A particle 
hops out of a randomly selected site i with occupancy 
Ui with a specified rate u(ni), and goes to site {i + 1). 
In the thermodynamic limit TV— >oo,L— >ooata fixed 
density p — N/L, the probability P{C) of a configuration 
C = {ui} in the steady state, in the grand canonical 
ensemble, is given by [3| 



integrated current H{i,t) which counts the total number 
of particles crossing the bond + 1) in time t. Thus, 
we have 
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if n > 0, 
if n = 0. 
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Here, Z is the partition function and v is the fugacity. 
The average mass at a site equals vF' (v)/ F{v), where 
F{v) = 1-1- ^'^^iv"f{n). Particle conservation gives a 
relation between p and v, namely, p = vF' (y)/ F(v). The 
fugacity has the maximum value Vmax = m(oo), given by 
the radius of convergence of the infinite series F{v). The 
ZRP can be mapped to a generalization of the asymmet- 
ric simple exclusion process (ASEP) by interpreting the 
ZRP sites as particles in the ASEP, while the particles 
at a ZRP site become holes preceding the corresponding 
ASEP particle [1]. The hop rate u(n) for an ASEP par- 
ticle, now a function of the headway to the next particle, 
induces a long-ranged particle hopping. 

Here, we consider the hop rate u{n) = 1 + b/n with b > 
2 for which the system undergoes a nonequilibrium phase 
transition As p crosses the critical value pc = 1/(6—2) 
0, a low-density disordered phase with mass of 0(1) at 
each site evolves to a high-density condensed phase where 
a macroscopic collection of particles of average mass {p — 
Pc)L condenses onto a randomly selected site, while the 
remaining sites have the average mass pc- 

In the steady state, the mean current of particles be- 
tween every pair of neighboring sites is the same. In 
the thermodynamic limit, the mean current is J = 
'^^^iu{n)v^f{n)/F{v) — v. It increases with the den- 
sity p, attaining its maximum value Umax = u(oo) = 1 
at Pc, and remaining pinned to this value in the con- 
densed phase. To address the dynamics of density fiuc- 
tuations, we examine the fluctuations in the integrated 
particle current across any bond in the steady state. The 
large deviation function of the integrated current has 
been studied for the ZRP with open boundaries in 
For our purpose, we monitor the variance C{i,t) of the 
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In the equivalent ASEP, C{i,t) is a measure of tagged 
particle correlations, being given by the variance of the 
ith tagged particle around its average displacement in 
time t. In simulations, we monitor C{t) = YlLi 0/-^ 
and find that it shows strong differences in behavior in 
the disordered and the condensed phases, refiecting very 
different underlying physical processes in the two phases. 

Disordered phase. C (t) in this phase behaves similarly 
to the tagged particle correlation in the ordinary ASEP, 
studied recently As discussed below, there are two 
scales Ti ^ L, set by the circulation time of a kinematic 
wave of density fiuctuations, and T2 '~ L^^^, given by the 
time taken by this wave to decay. 

(i) t <^Ti. Here, C{t) grows linearly in time: C{t) — 
vt. This follows from the result that, in this time regime, 
H{i,t) is Poisson distributed with intensity v over all 
bonds 1). The population at a ZRP site under- 
goes a time-reversible birth-death process where a birth 
(particle input) occurs with rate v, while a population of 
n particles undergoes a death with rate u{n). For a re- 
versible birth-death process with Poisson inputs, Burke's 
theorem implies an identical Poisson distribution of out- 
puts [13] ■ Noting that the output from one site forms the 
input to the next site then implies the result [11|. 

(ii) Ti <^ t T2. In this regime, C{t) oscillates as a 
function of time. In a driven system with homogeneous 
density p and a density-dependent current J(p), density 
fiuctuations are transported as a kinematic wave with 
speed VK = dJ/dp [l2|. This wave is dissipated over a 
time scale ^ , where z is the dynamic exponent of the 
system. For the ZRP, z takes on the Kardar-Parisi-Zhang 
(KPZ) value of 3/2 [l3|. Since z > 1, fiuctuations circu- 
late several times around a periodic system before getting 
dissipated, and revisit every site after a time L/vk- This 
makes the variance oscillate in time with this period. A 
measure of the growth of dissipation in time is given by 
the lower envelope of the oscillations, which behaves as 
^2/3 where /3 = /3kpz = 1/3 [H. 

(iii) t ^ T2. The time scale T2 marks the dissipa- 
tion time of an initial density profile. For times t ^ T2, 
the variance grows diffusively: C{t) ^ D{L)t. Match- 
ing this behavior at T2 with that in (ii) above gives 
D{L) - as for the ordinary ASEP 

Critical point. At the critical density pc, the variance 
behaves differently for values of & < 3 and & > 3. For 
6 < 3, there is no moving kinematic wave 0]. Hence, the 
integrated current is Poisson distributed with intensity 
'^max = 1, implying that the variance continues to grow 
linearly with slope 1. For 6 > 3, however, the kinematic 
wave speed is nonzero and the Poisson distribution for the 
integrated current is expected to hold for times smaller 
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FIG. 2: (Color online) Monte Carlo (MC) simulation results 
for B{t)/t^^^ as a function of time t for different fe's at pc — 
l/(b — 2). The system size is L = 16 384. The dashed lines 
show the asymptotic t^^^ approach. 
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FIG. 3: (Color online) (a) Schematic plot of C{t)/t as a func- 
tion of time t in the condensed phase. The various regimes 
(see text) are also marked, (b) Scaling of C{t)/t with t/L'' for 
6 = 3 in the condensed phase (p = 4). Various system sizes 
L considered are marked in the figure. The data points, ob- 
tained from MC simulations, are connected by smooth curves. 



than the return time of the kinematic wave. At critical- 
ity, the largest mass in the system scales as ^ L'^/i^-i) 
[ISj and is insufficient to block the circulation of the kine- 
matic wave around the system. C{t) oscillates in time as 
for p < Pc with return time Ti and decay time T2 of the 
kinematic wave. To find the exponent /3 at criticality, we 
monitored the variance B{t) of the integrated current by 
starting from an arbitrary but fixed initial configuration, 
drawn from the stationary ensemble [l6| . B(t) grows 
asymptotically as t^/? for t < 2.3/2 ^ f^^^ ^^lat /3 
at criticality has the KPZ value of 1/3, independent of b 

(Fig. m- 

Condensed phase. For p > pc, a finite fraction of the 
total mass (the condensate) resides on one site for the 



characteristic survival time Tg ^ {p — p^^^^L^ [17|. It 
then relocates to another site over the relocation time 
scale Tr ~ {p — pcY'L'^ ^ as discussed below. The behavior 
of C (i) in this phase is best depicted by plotting C (t) /t as 
a function of time, as shown schematically in Fig. [3fa), 
where the various regimes are also marked. 

(i) t <^Ts: Here, C{t)/t equals 1, with a mild upward 
deviation for longer times, (ii) t Tg: In this regime, 
C(t)/t rises rapidly in time, (in) t>Ts+Tr: Here, C{t)/t 
falls slowly in time, (iv) t ^ Tg+Tr'. Here, C{t)/t begins 



FIG. 4: (Color online) (a) Survival probability distribution 
Psij) of the largest mass in the condensed phase. The sys- 
tem size L = 128, the parameter 6 = 3, while p = 4. The 
dashed line is a guide to the eye for the part of Psij) behav- 
ing as T~^^'^ . (b) Scaling of the bump in Ps{t). The data are 
obtained from MC simulations. 



to approach a size-dependent constant, as predicted by a 
simple model described later in this paper. Features (ii), 
(iii), and (iv) result from enhanced fluctuations due to 
the condensate relocation. To understand this, we need 
to first discuss the relocation dynamics. 

The condensate relocation occurs through exchange of 
particles between two sites. On monitoring the time evo- 
lution of the largest and the second largest mass in sim- 
ulation, the following picture emerges. Let M{t) denote 
the largest mass in the steady state at time t. M(t) has 
the average value Mq = {M{t)) = {p — pc)L, with fiuctua- 
tions AMq which scale as L^/^ for 6 > 3, and as L^/(''~^') 
for 2 < 6 < 3 [18|. These fluctuations may build up in 
and over the time scale Tg, the largest mass fluc- 
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tuates to ~ Mo/2, while a mass ~ Mo/2 also builds up 
at another site. Subsequent to this, two sites with mass 
^ Mq /2 exchange particles between themselves resulting 
in relatively rapid alternating relocations of the largest 
mass from one site to the other. The difference of masses 
on these two sites performs an unbiased random walk in 
time until fluctuations populate one of the sites to ~ Mq 
at the expense of the other, which happens over the time 
scale Tr {p — pc^L'^. 

Figure [4l^a) shows the Monte Carlo results for the sur- 
vival probability distribution Pg{T) of the largest mass, 
obtained by computing the distribution of the time inter- 
val between successive relocations. Pg{T) has two parts, 
(i) a power law part ^ r~'^/^ and (ii) another part, which 
corresponds to the bump in Fig. HJ^a), and has the scal- 
ing form (/9-pc)~*''+^^i"^^+^V i-r/Tg), as shown in Fig. 
m^b). The prefactor comes from the normalization of 
Ps(t) to unity, with the cutoff for the t"^/^ part taken to 
scale as T^. The power law part holds for times when the 
two sites with mass ^ Mo/2 compete to hold the largest 
mass. The random walk argument of the preceding para- 
graph predicts a t"^/^ decay, since Ps(t) then stands for 
the probability for the random walker to cross the ori- 
gin for the first time. The second part in Pgir) arises 
from the relatively long time for which the condensate is 



stationary on one site. 

We now explain the behavior of C (t) /t in the different 
regimes seen in Fig. |3lja). The condensate is stationary 
on one site for a long time ri , which is a random interval 
distributed as p{n) ~ {p - pc)-^''+^'> L-^'>+'^^ f (n/Ts), 
with the characteristic survival time T,. In regime (i), 
when t ^ Ts, the condensate is stationary and acts as 
a reservoir for fluctuations, preventing their transport 
around the system as a kinematic wave. As a result, 
Burke's theorem is valid over this time scale and the 
integrated current is Poisson distributed with intensity 
I'max = 1, implying C{t)/t = 1. When t ^ Tg, the con- 
densate starts to move from one site to another by trans- 
ferring its mass across the intervening bonds. As a result, 
these bonds pick up enhanced fluctuations [oc (AMq)^] 
in the integrated current over the relocation time inter- 
val T2 , which is a random variable with the characteristic 
time Tr- These enhanced fluctuations lead to the rise 
in C{t)/t as a function of t in regime (ii). The collapse 
of the rise times seen in the scaling plot of Fig. [3)^b) 
confirms this picture. In regime (iii), after t + Tr, 
the condensate has completed relocating, so current fluc- 
tuations revert to Burke-like behavior, resulting in the 
fall of C{t)/t in time. The slow fall in regime (iii) arises 
from the wide distribution of the time T2 , and further re- 
locations. To predict the behavior in regime (iv), where 
t 3> Tg+Tr, we construct below a simple relocation model 
(RM) that describes the effect of condensate relocation 
on the long-time behavior of current fluctuations. 

Figure [5] shows schematically the instantaneous cur- 
rent t) across the bond (i, i -I- 1) as a function of time. 
In a given history, let K be the number of condensate re- 
locations in a fixed time t. Here, if is a random variable 
with mean given approximately by {K) « t/{DTs+BTr), 
where B and D are constants, independent of the density 
and the system size. At the kth relocation of the con- 
densate (fe = 1, 2, ... , K), let Mk denote the amount of 
mass transferred across the bond (i, z -I- 1) over the inter- 
val T2. The random variable Mk has an average Mq and 
variance AMq. The integrated current H{i, t) is given ap- 
proximately as H{i, t) w Y^k=i + /g dt'j{i, t'). 
For & > 3, on computing the variance of the integrated 
current, we get C{t) w GL{K) + [t - (K)BTr), where 
G is a constant. On substituting for (K) and neglect- 
ing the time scale Tr ^ (p — PcYL"^ in comparison to 
Ts ^ [p — pcf'^^L^ , we obtain the asymptotic behavior 



C{t)^[L-%p-p,)-'^^+^^ + 1% 
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with 9 = b — 1 for 6 > 3. For b between 2 and 3, we find 
that 9 = {b'^ - b- 2)/{b - 1). Thus, the RM predicts 
that, for all values of 6, at long times t ^ Tg + Tr [regime 
(iv)], C(t)/t approaches a size-dependent constant which 
scales down with the system size. This long-time regime 
(iv) could not be accessed in simulations for the system 
sizes shown in Fig. [Sfb), but we confirmed its existence 
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FIG. 5: (Color online) Schematic plot of the instantaneous 
current j{i,t) across the bond + 1) at time t as a func- 
tion of time. The random variable ri is the time for which 
the condensate is stationary on a site, while the random vari- 
able T2 stands for the relocation time of the condensate. Mk 
is the integrated current over time T2, arising from the fcth 
relocation of the condensate across the bond + 1). 



for smaller systems. 

In summary, we addressed the dynamics of steady state 
fluctuations of a zero-range process which undergoes a 
nonequilibrium phase transition from a disordered to a 
condensed phase. Different dynamical properties emerge 
in the two phases. In the disordered phase, fluctuations 
move around the system as a kinematic wave. Such a 
wave, though present in the bulk, cannot circulate around 
in the condensed phase because the condensate subsumes 
fluctuations. The dynamics is governed by the conden- 
sate relocation through a slower process of transfer of 
particles from site to site, contributing enhanced fluctua- 
tions to the particle current across the intervening bonds. 
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